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ABSTRACT 


NASA has recently published data on the effects of ten- 
perature and pressure distortion on the stability of axial- 
flow compressors. Theoretical predictions of compressor 
instability based on vorticity do not explain the observed 
compressor stability of the NASA tests. Description of 
distorted flow by vorticity is extended by a simplified model 
which incorporates vortex sheets and influence coefficients. 
To determine the upstream influence of the compressor on 
flow field distortion, a,model using parallel compressor 
theory, experimental data from NASA and a pressure ratio 
equation derived by Hill and Peterson is formulated. The 
model correctly predicts compressor blade loading calculated 
by Carta. This new model accounts for upstream effects of 
the compressor and shows that the vorticity model by itself 
is incomplete. Incorporation of the new model into a dis- 
tortion index based on radial vorticity and total pressure 


change is proposed and discussed. 
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I. INTRODUCTION 


In applications of high-specific-energy turbomachines, 
such as turbo-jet engines with axial-flow compressors, it 
is necessary to operate the compressor as close to maximum 
pressure ratio as possible without inducing unstable per- 
formance known as surge. A phenomenon of axial flow con- 
pressors is sensitivity to non-uniform or distorted flow 
fields when operating near maximum pressure ratio. As air- 
craft are designed for higher speeds and more comprehensive 
operating envelopes, the possibility of flow field distor- 
tion at the compressor face is increased. At supersonic 
speeds, shock effects within the engine intake system are 
significant. At high altitudes and low Mach numbers the 
effects of low Reynolds number may cause portions of the 
flow in the compressor to revert to laminar flow.  Laminar 
flow increases the tendency toward separation which accel- 
erates compressor stall. Aircraft operations at high angles 
of attack and/or yaw also can cause distortions of the flow 
7701 nip performance military aircraft may be subjected 
to all of the above distortions in addition to armament re- 
lated temperature distortions. Forward firing armament 
could also cause ingestion of foreign gas into the compressor 
and is discussed in Appendix A. Compressor surge is not re- 
stricted to the airborne flight envelope. Ground operations 


with some recently developed engines have resulted in 


iL 





destructive surge while taxiing. Fıgure 1 shows some of the 
general areas of the flight envelope where large stability 
margins are required. Figure 2 shows a compressor map that 
illustrates the stability margin. The associated performance - 
penalty of a large stability margin becomes more acute as 
engine performance requirements increase. 

This situation has motivated the formulation of a dis- 
tortion index that could be used in the design, development, 
and component matching of turbo-jet engines which use axial- 
flow compressors.  Conceivably such an index could be used 
in engine control systems to allow operation at reduced 
stability margins without the danger of compressor stall. 

Several approaches to distortion indices have been tried. 
In fact there has been a proliferation of indices based on 
pressure distortion or temperature distortion at the face 
of Lhe compressor. Most data collected are in the form of 
total pressure or total temperature measured at the com- 
pressor face. A more recent approach is to formulate an 
index based on vorticity. Since vorticity can be derived 
from either temperature or pressure data, a single index 
appropriate for both forms of distortion is possible. 

Farmer [Ref. 1] investigated converting total-pressure 
Feen: ars the compressor Taceı Into vortieity. Iverson 
[Ref. 2] proposed a method of transforming total-temperature 
5 if ürS into vorticity.” Shoemaker [Ref. 3] expanded Farmer's 
theory and, using a computer generated map of vorticity, de- 


Eu dan index, based on vorticity alone, that predicted 
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stall in selected stall events. Research conducted, concur- 
rent with this paper, by Vann [Ref. 4], refines the Shoe- 
maker vorticity index and applies it to additional stall 
events, The refined vorticity index (Vann) is an accurate 
Stall predictor in 45 percent of the cases studied. 

Braithwaite, Graber and Mehalic [Ref. 5] address the 
Eme cts ote temperature and pressure distortion on 
turbo-jet performance and provide experimental results with 
which the vorticity model can be compared. In the case of 
boralenressure dıstertion, tne vorticity model is a good 
predictor Of Gcereased Compressor performance. INe= vert Lele 
model based on total-temperature distortion does not cor- 
relate well in all cases when compared with Ref. 5. 

Toca further insight into this-apparent failure tor 
Ghe vorecaty model to correlate with experimental results, 
this paper will apply the vorticity model to models sımu- 
lating the conditions of Ref. 5. "SEnwaddrtronswNancextenscis 
model that includes the engine effects on the upstream flow 
field will be introduced and a new approach to a distortion 


index will be discussed, 
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EEE NATURE OF THE PROBLEM 


The vorticity approach of Farmer, Iverson and Fuhs 
[Ref. 6] gives the following relationships for pressure and 


temperature distortion: 


] aP ' 
u) 1 MN ME Erro m “ERRE (1) 
T U YP T 90 
M Su 
Or ^ 2r" 305 = 


where variables with superscripts (primed) are non-dimen- 


Sionalized quantities as defined below: 


r = r's D = compressor diameter 

UE a - average speed of sound 

wW = ,'28 w = fluid vorticity 

ри an = average stagnation pressure 
T. = T I - average stagnation temper- 


ature., 


иси ир T r denote circumferential and radial 
1405000016  Tespectively. A complete list of symbols is in- 
sır ci ompages 8 and 9. 

571600  ationsBbıps given in Ref. 6, it also can be 


shovn that in the case of combined pressure and temperature 





distortions 
3P 1 19T 
7:7 Uy .. ta. E с un (3) 
LU rz гР, 290^ iA ‘Gan! M? pro : 
2 





Order of magnitude analysis [Ref. 6] indicates that Uw, is 
of negligible size. Reference 5 compares the effects of 
m sure and temperature distortion, separately and in com- 
bination, on the performance of J85-GE-13 turbo-jet engines. 
To induce total-pressure distortions on the face of the 
compressor, screens of varying solidity and 180” circumfer- 
ential extent (semi-circle) vere placed one diameter upstream 
from the compressor. The total temperature of the air vas 
distorted by means of a hydrogen burner, also of 180” extent, 
located four diameters upstream from the engine inlet. The 
burners vere external to the bellmouth. The test procedure 
is discussed ın detail in Ref. 5. Figure 3 shows a simpli- 
fred schematic of the test set-up used. Figure 4 is a com- 
pressor map showing the degradation of performance associated 
with three combinations of pressure and temperature distor- 
tion at the engine inlet. From Fig. 4 the most adverse 
effect on compressor performance is observed to be the case 
of temperature and pressure distortions super-imposed (180° 
overlapped). Figure 5 shows the circumferential profiles 
Dto tel temperature and total pressure at the engine inlet. 
Based on the vorticity approach to distorted flow, one 
would expect that large values of oM would correlate with 
re loss of stall margin. Radial vorticity can be calcu- 
İB (roy Edn. (5). When 9P,/30% and 9T. /30* have the same 
algebraic sign, one would expect the most unfavorable case. 
Reference to Fig. 5(a) shows that for the 180° overlapped 


| case, 3P, /30* and 9T, /36* Mavesoppositessnens:; Tor thisMcase 


1:7 





“J9TUT 1osseiduo) IY PISTA MOTA UOTINOIST( asoduj 03 pssn dn-39S 259] Fo AtJeusyas  'c oin$ijJ 





dv 
— – 
q 


u3nouTTeq 
Iossoiduo иәәлос̧ Iouing 








'(S "yeu uoig) UO0TI1OIST( sın? 
-eleduo[ pue oinssoig peurquo) 3uoiXxg 4,081 jO S3293JjH SurMoug dey 1iossouduo; 


2 | 
INOA “(o/eA®m) 'mo1salv 021938809 
201 86 16 06 98 28 8/ m 0) 





IN2IE eo 
Coh hN \ ¢ 


(3345 10103 og 
0312933802 Y 


001 
| 
\ 


NOWROLSIUBON —-— 


(344V133A0 g081 (04 v 
1344Y11340 (06 (24 a 
> 0IAAVIYVIAO (0 18d О 


a | NOILYOLSIG 





“y Janets 


CdIEd 'OLLVY JYNSSIUA YOSSIUAWOD 


19 





O TEMPERATURE 
1.10 о PRESSURE 
„--{P=---- l ke =, 







= 1.00 

<L 

oa 

Zw 90 

EI (A) DISTORTIONS 1809 OVERLAPPED. 
— 

ui 1.10 ! 

= x—-----0----. 

E 

ə 

ud 

oa 

= 

”“ 

UN 


INLET PRE 


0 90 180 270 360 
CIRCUMFERENTIAL LOCATION, 8”, DEG 


(C) DISTORTIONS 0% OVERLAPPED. 


Figure 5. Compressor Inlet Temperature and Pressure Profiles 
with Combined 180? Extent Pressure and Temperature 
Dxstortions (From Ref. 5). 
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o. is small, and there should be only small loss of stall 
Margin. For the 0° overlapped case, Fig. 5(c), 9P,/90* and 
9T, /30* have the same sign, and W is large. The 0° over- 
lapped case should result in large loss of stall margin. 
Observing Fig. 4, it can be seen that the opposite result 


was observed experimentally. 


21 





III. MODELS OF FLOW WITH PRESSURE AND TEMPERATURE DISTORTION 


Bi Tects on tir flow field at the face of the compressor 
when distortions are imposed upstream by screens (pressure) 
a tano (temperature) will be examined. Tables contains 
a list of influence coefficients for area change (dA/A), heat 
addition ED and #riction (dC £/2) that will be used in 
this analysis. These influence coefficients are adapted 


epom ohapiro [Ref. 7]. 


A. PRESSURE DISTORTION 

Figure 6 15 a schematic of the modell usedinon pressure 
distortion. The streamline curvature depicted is for 
graphic purposes, and it is not assumed that the streamtube 
eMange actually will be Ms @iawn. “The sereen is assumed to 
be of 180° extent (which is consistent with the methods of 
Leo) sancdethe flow at Station 1 is assumed to be uniform: 


From continuity 


Maz - m = pU A = mass flow rate 


A2 I 


2:53 37 £eTeN64al form düter TOZyarithmic dıdferentiation 
"= ј= 970, 40, ЈА 
dm = 0 = 5 E ae 
The physical geometry requires that 
dA dA 
Da "Wa: 


The assumed uniform static pressure distribution requires 
that 


po 
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AREA CHANGE 





dA — 
A C ии 
p 

dU 1 - 1 
U M^ - 1 M2 - 1 
B Y. N yM? 
p M? - 1 M? - 1 
E . 
T Y - 
t ın... 
TABLE 1. 


ander 1 Ct lon : 


dP, _ dB, 
P’A DB 
From Table 1 

qıs =, sels dA, 
P’B M2 - 1 А”В 


dP, ME 
P-A M2 - 1 


HEAT ADDITION 





dA, 
AA 


FRICTION 
d Ce 


2 


_ [1 + M^ (y - 1)]y M? 
Tam om 


Influence coefficients for area change, heat addition 


(5) 
(6) 
(1 * М2(у-1))у М? Tr 
jM ZN: 
(7) 


Combining Equations (4), (5), (6), and (7) and algebraic 
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manipulation yields the following relationship for area 


change: 


Ал = (14 М2 (у - 1)1 = - 90. (8) 


Wes Similar techniques to investigate velocity, it®can be 


shown that 
dU, . — 1. dA , yM*' (Tr . 
U A M2 - 1 AGE 1 - M? 2 | 


and combining Equations (8) and (9) with some simplification 


reduces to 


dC 
ducc Fiat £ 
WA | 2(1 - Mi | pee om 
From Table 1 
dU, 1 dA 
NR 1-N a (11) 


and the ratio of velocity changes can be expressed as 


al, Beam 1 dog 
A 2 2 
gU. o E EAS 

U- B E 11- A'B 


which, using Equation (8), reduces to 


dU 
mA LIMINE 12:4. M* = 1 . 
Me = ==. + 50) dz 


U 
“Әв 


for y = 1.4. 
Equation (8) shows that for all real valued Mach numbers, 
Streamtube A will increase in cross sectional area (the fric- 


tion coefficient is always positive), and streamtube B must 
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decrease a like amount as in Equation (4). For subsonic 
flow, Un 2 must then be greater than Uni: Equation (12) is 
finite for all real values of Mach number but changes sign 
when M = 0.6455. For values of M less than 0.6455, Equa- 
Egon (IZ) is negative which@requires the velocity Ghange in 
streamtube A to be negative. Conversely when M > 0.6455, 

U 


> U Quantitatively M may be less than 0.6455 al- 


A2 A1’ 
though transonic compressors have axial Mach numbers ap- 


proaching 0.65. 


B. TOTAL TEMPERATURE DISTORTION 

Figure 7 is assehemetic of thewtotal temperature déstor- 
tion model. The illwstrated curvature”of the streamtube is 
for graphic purposes and does not represent an assumed pat- 
tern of change. At station 2 the static pressure distribu- 
tion is assumed to be uniform. This implies that station 2 
1s far enough ahead of the compressor so that there is no 
HE canub or the engine. 


Using Table 1 


dP . 2 YmM dA 

HU) p = l - M2 A)B e) 

dPj . yM* dA) . yu Ün (14) 
P-A :2—- AA 


1:51 ии 
P 


Setting Equation (13) equal to Equation (14), because of the 


uniform pressure assumption, then 
dh 


dA ОЛ if 
g | > pT (15) 
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Since 


dA, = - da, 
A B AA 


Equation (15) becomes 


dh 
dA dA 1 E 
P 


Wasch the addition of heat, SN decreases and SA increases; 
and for subsonic flow, Un) must be greater than UR: 

The effect of heat addition on the velocity is investi- 
gated similarly 


1 dh 
IP CT (17) 
p 


m. 


dU, 01 dA 
"MEN ил DA? 


dU dA i 
mmm M” US 


Incorporating the streamtube area changes of Equation (16) into 
Equations (17) and (18), the velocity change is found to be 
the same for both streamtubes. 

1 dh 


dU. x dü, x^ gee 
UA^ UB" 2(0-H cT 


(19) 


C. COMBINED TEMPERATURE AND PRESSURE DISTORTION 

Figure 8 is a schematic of the models used for combined 
pressure and temperature distortion. Cases will be considered 
for zero degrees overlap and 180" overlap. 

From Table 1 the addition of heat causes a loss of stag- 


nation pressure. 





4Р, ey M? /2 dh, _ Y M? dh, 
| 5 1423 M? Es 2 A 





“dETİƏAQ 4,081 PU? 4,0 103 TOSPON UOTILOISTA DANSSIAG pue sınJeilsdus] poutquoy 


“dETIƏAQ 4,0 Sinsssig pue sınyelsdusL (q) 





d 
Iossoiduo-; UƏƏZƏS 
X 
x V 
o08I o0 ТК ра ee ee 
199H 
*dBTIDAD ¿08T SANSSIAG pue sinjeredusay (©) 
2 c 
10ss91duon U9S1IS 3?9H 


o08I 0 





2107000011 


27 





In the case of 0° overlap 











t - Y.M” 
md Г— ас̧ > 
S ə M? E (215 
P. ^B TU OWN 


The net area change of streamtube A is, using Equations (8) 
and (16), 
dC dh 
dA f 1 
-— ) ANET ә [1 m lo) M?) 4 i p i = (22) 
To obtain results for a specific case, let us assume that 


ehe heat amd screens are such that 





dC. A dh, 
2 2 Colt 
Consequently 
dP dP 
E. cT s 
t t 


Foauation (22) in terms of l and introducing dCe/2 in terms 


of heat becomes 











JA, - [1 (у-1) M?) İİ TA 
AT ANET 2 > © T 
pr 
dh 
hee n 
pt 
A “ANET 4 uc 2 2 


veloc Change in streamtube A, 
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) = 1 | dA ) + m... M* -— 
U -A 1 - M? A “ANET e 2 


and incorporating Equation (23) simplifies to 


Hj a d d = (24) 
UB (1-M)j2 2 | 
In streamtube B, 
MB 1 - M2 A “BNET l- M? Cp! 
II 1- M? AAN 1 - M? En 
dC 
dU 3/2 + - 1)M? i 
mu ет 1 o : X c (25) 


Иһеттатлодог velocity hencges in the two streamtubes (for 





y = 1.4) xs 
zu s 
E E o M 
ume 05 a N) (26a) 
BU Y 


p -mudtaoneotobauatron (208) shows that the velocity in 
5177 0110 76 3 CBC for all subsonic values of M. Recall 
Buchen spes for the specific case of equal total 
2 T 557 nE 7 Dua Cİ resulted from the ar- 
570617 “0716 don of dCe/2 - ——.. 

For the case of 180” overlap the same method of develop- 
uni 0  Eigure 8(а) 15 the model used, and there are 


5:08 001306 illustrated for the purpose of adding the 


Sons ot each distortion. Assume that static pressure 
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is uniform at stations 2.and 3. Adding Equations (19) and 
(10) and again considering the specific case of 


dh 


nt 


ct 


dC, 
2 





- 1 
2 


the velocity change in streamtube A is 


a“ Ју = döz (1 ә 6) 
4(1-M2) 


Fer streamtube B, using Table 1, Equations (8), (16), and 


(19) the velocity change is 


dC 
dU E 2 
"uy Jp z 4(1-M2) e+ ae) 


and the ratio of velocity changes is, for y = 1.4, 


dU 


U IA _ 1+ 2 M? _ 1l + 2.8M2 _ 
dU, a a ” o (26b) 
U ^B 


Examination of Equation (26b) shows that the velocity in 
streamtube B is always higher than the velocity of stream- 
ube A. For the case of IS] aaa o distortions, the 
total pressure will not remain equal in both streamtubes 

as in the case of 0° overlap. In the case of 180° overlap, 
the change in total pressure for the distorted streamtube 
(A) will be 

dC e 


— a ә 


2 


-2YM? 


for this specific case of 








Total pressure will remain constant in the undistorted stream- 


tube (B). 


D. HEAT ADDITION AT SMALL MACH NUMBERS 

Using the temperature distortion model of Figure 9, the 
effects of heat addition on the flow field will be made more 
mestétrictive by limiting the entering Mach number to very 
small (near zero) values. This is accomplished by placing 
the burner outside the bellmouth, where velocities are small. 
This analysis also models the tests by NASA reported in 
Ref. 5. An energy balance at stations 1 and 2, assuming 


adiabatic conditions between these stations, gives 


I 12 UE = 
Ne) ти (ии ил 


; l y? i l y2 
17 e U 22 * mp Be + 5 UR) CC. | 557 


Note that station 1 is immediately dovnstream from the burner. 


For streamtubes A and B then 


UNES 02 270 DA 57 
and 

75. ve = cem e DP. secum (29) 

P BI 2551 p B2 zer? PoBo : 
Assume uniform P at stations 1 and 2, Ma = mp and sent opus 


flow for streamtubes A and B between stations 1 and 2. With 
the previous assumption of small Mach number, conditions at 
station 1 are specified. The stagnation pressure for both 


Are amtabes at ertier Station 1 or 2 equals the ambicnt 
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Figure 9. Temperature Distortion Model for Heat Addition at 


Low Mach Number. 
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pressure. Static pressure for both streamtubes is equal at 
station 1. Due to heat addition Tra > Tip: The following 


isentropic relationships now can be used: 





5 Nasa | 
575-— E. | — 
and 
| у:1 
152 7 Teg p— PE | (31) 


Using Equations (30), (31) with (28) and (29) will give the 


relationship of velocity for streamtubes A and B 


Y-1 
p Y 
U = xT 1-12 
A2 p tA P. 
у-1 
Р 35855 
ne cT 1 - 1-2 
B2 p P. 
2 
“255. (32) 
U ik 
B2 CB 


Since 


P 
oJ, = pops (1) 


and for hydrogen (fuel for burner in NASA test of Ref. 5) 
i 01296 in exXxcess air, Ya 7 Ypg» then Ma = My - 5755 
Investigating the area change at 2 with the same assump- 


tions 
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pAU),, = pAU),, => Gee = pe pt 3 
B2 A2 B2 VT 
tA 
2 2 
5050, 
o B2 AD 
and using Equation (32) 
10, | Ne 
A 
B2 YT 


Equation (34) is restrictive in that the assumptions limit it 
to the case of equal mass flow in streamtubes A and B. Hov- 
ever, parallel compressor theory uses a mass-weighted sum of 
subcompressors to determine the operating point of a distorted 
compressor, so in the case of 180? sector distortion this as- 


sumption is not unreasonable. 
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IV. COMPARISON OF FLOW MODELS WITH VORTICITY THEORY FOR TEM- 
PERATURE DISTORTION 
The analysis of Section III(B), total temperature dis- 
tortion, shows that the velocities of the two streamtubes 
are the same, thus making the radial vorticity zero. Set- 
ting o in Equation (5) equal to zero and solving for oP, / 30% 


yields the following relationship: 








T 1 
x * = 
TP.U, 98 rU, 9 0 1 4 Y M? 


RT, Py Yl oae 
P. p Lo 
dP _, dh c 
yocant 
t i p 








S 


5111177077) s equivalent to the influence coefficient of 
1777077) 11005 verdiying VOTtiCity Гот this case. 

For the case of heat addition at low Mach number, a more 
subjective comparison must be made.- Equation (3) simplificd 


for temperature distortion only shows 


5 





(36) 


x 


ct 
— 
1j 062 


o 
Qo 
=] 
ә 
1 
EA 
a a 
N 


S 


— 
T RU, 99* 


All values in the right hand side of Equation (36) are always 
positive with the exception of the circumferential total tem- 
perature gradient. For the velocity profile shown in Figure 
9, w, is positive when 0* = 180° and '9T,/30* is valisso posi- 
tive. A similar comparison at 6* = 0° also shows algebraic 
agreement in that 0). and the total temperature gradient are 


both negative. 
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V. EXTENDED MODEL TO ACCOUNT FOR COMPRESSOR INFLUENCE 


A. PARALLEL COMPRESSOR MODEL 

Reference 5 employs the parallel compressor model as a 
method of predicting the effects of distortion on stall 
pressure ratios. Figure 10 shows a schematic of the parallel 
compressor model with some of the assumptions used. Figure 
11 shows the predicted stall points of the compressor for 
various pressure and temperature distortions using this 
modell Figures 12 and 13 show the profiles of static pres- 
Eure, total pressürerand total temperature at the compressor 
inlet and outlets for both pressure and temperature distor- 
tion. From Figure 12 it can be seen that static pressure 
is not uniform at the face (inlet) of the compressor. This 
was one of the assumptions used in the development of Sec- 
mom lit or thas paper. The development of an extended 
model in this section will consider the effects of non- 


uniform pressure profiles. 


E PRESSURE RATIO EQUATION 
Hill and Peterson [Ref. 8] derive a compressor pressure 
5707 -:2 :030n that predicts the total pressure ratio (N) 


across an axial compressor stage as 


2751: 0) nu. — e 
E | SE tan 81 


(37) 


where the subscripts 1 and 2 refer to stations before and 
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|. SECTORS C AND D ARE CONSIDERED AS 
INDEPENDENT PARALLEL COMPRESSORS 

2. BOTH SECTORS FOLLOW SAME MAP AS 
UNDISTORTED COMPRESSOR 

3. NO CROSSFLOW BETWEEN PARALLEL COM- 
PRESSORS 

4 BOTH EXIT TO A COMMON STATIC PRESSURE 


erect Parallel Compressor Model (From Ref. 5) 
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Figure 11. Application of Parallel Compressor Model to Predict 
: 757511 Fo:ints (Prom Ref. 5). 
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Figure 12. Compressor Temperature and Pressure Profiles with 
180" Circumferential Temperature Distortion (From 
Ref. 5). 
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STATIC PRESSURE, TOTAL PRESSURE 
OR TEMPERATURE RATIO 
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(B) COMPRESSOR EXIT. 


Compressor Temperature and Pressure Profiles with 
180° Circumferential Temperature Distortion (From 
REUS. 
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after the compressor stage. Figure 14 shows the relationship 
of the angles By and. В. with a compressor blade row at wheel 
speed Qr. For flow in unstalled cascades of high solidity 
(ratio of blade chord length to blade spacing in the cascade), 
the exit angle B, is relatively constant with changes of Ву. 
Blade loading is a funetion of, or can be correlated to, AB. 
If -B, ~Gons tant , «then B4 has a significant effect on the 
stall characteristics of a cascade. Consider Equation (37) 
with a fixed n and wheel speed. An increase in temperature 
must be accompanied by an increase in B1. An increase in B1 
can be obtained only by a decrease in axial velocity. This 
is shown in Figure oe Figure 11 shows that parallel com- 
pressor theory substantiates the decrease in axial velocity 
for this case since corrected weight flow is proportional to 
axial velocity. “The model of Section IJI.D. shows a veloc- 

yı “012127” TeOV€TSsCÜ Iron the Tesulgszor this Section; see 
Figure 9. An assumption of Section III.D. was uniform static 
5257551055) the downsticam station. Ihe model in fact pre- 
cluded any upstream effect of the compressor itself on the 


flow field. 


C PRESSURE DISTORTTON 

1710110700 l lustrat es chesprosiles of static and total 
pressure, Static and total temperature and axial velocity 
as derived from the models of Section III and as measured 
at the face of the compressor during the NASA tests Oa cue 


ENMENDIGUVEDPODpOrtres atothe constant-static-pressure station 
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agree with those illustrated by Figure 6. According to 
Figure 12, the static pressure increases in the high-P, 
streamtube and decreases in the low-P, streamtube. The 

flow decelerates in the Һати streamtube and vice-versa 
for low-P, streamtube. Also the area of the high-P, stream- 
tube of Figure 6 must increase. 

Me screen creates a vorticity sheet which is convected 
into the compressor. This sheet separates the high-P, and 
low-P, streamtubes. As predicted by Helmholtz's vorticity 
Wows the vortex sheet 15 a flow feature which persists once 
it is created [Ref. 9]. The vortex sheet decays by means 
of viscous forces. Figure 17 shows a model that uses a 
"bent vortex sheet" associated with the change of velocity 
mid streamtube area for pressure distortion. Specific points 
2575710) in relation to this vortex sheet are indicated, 
and the change in axial velocity at these points is also 
shown in Figure 17. 

DUCUM OG pende or tner voten meet 1167 lon acquires a 
Ug component in the direction of rotor motion at points A 
and B and a Un Opposed to TIGI R C and D of 


Figure 17. This has an influence on 81 as discussed later. 


D. TEMPERATURE DISTORTION 

Ss Mos ic £lomsproperties at three stations: 
immediately downstream of the burner, at a constant-static- 
pressure location within the duct and at the rotor inlet 


Station. The flow properties at the constant-static-pressure 


47 





|| | 


05 





Low PL 


(Screen) 


B! 


Bent Vortex Sheet 


AU_ Due to Area Change 


2 
D d A D 


AU, (*) | 
E 
0 = 
E 
0 90 0 
= 
B 
C 180? İ C 
AU, Due to Vortex Sheet 
Bent Vortex Model for Pressure Distortion. 


48 





'UOI3J103SIQ ednssoJq OJI97 U3TM UuOI3J1O3SI(q 9oJn3?u1oduo] “gt sinsty 


jerug UOI39?3Sg oJunssoJg 
IO30M IV 02t3464077108751602) “ON UOEN MOT ALOA 
“aad on "L ^a n a "aa ın 
Ol m m UOT31PPV 
3*9H 
FA FOU 


Jg 


49 





station within the duct agree with Figure 9 of en Ill 
and figure 13(a). 

177001 € 10 Cum in iluence of the compressor, 
focus attention on Figure 11(b) vhich is the parallel com- 
pressor model matching Figure 18. The hot streamtube, 
identified as 1 in Figure 11(b), has sharply reduced cor- 
rected weight flow. This implies lov axial velocity com- 
pared vith the cold streamtube, number 2 in Figure 11(b). 
17 11600) hc velocity profile reverses from fast in 
the hot streamtube at the constant-static-pressure station 
102) l1lov in (he hot streamtube at the rotor. This reversal 
s attributable to upstream influence of the compressor. 

The original vorticity sheet generated by heat addition re- 
mams. the static pressure is not uniform at the rotor, 
17771 vaş an essential element of the models of Section 111, 
битә ти TİC compressor müst be to change the areas of 
the streamtubes by an increased corrected veight flov in 

the undistorted streamtube. Figure 19 shows the bent vor- 
tex sheet as applied to the temperature distortion model 

and the velocity changes at the rotor. 

175751 1 707 TEE BENT VORT X SHEET ON ROTOR BLADE INCIDENCE 

ANGLE 

итә -:mccy0 Rer TO] postulate that it is not 
бә ти reseribe the conditions at the compressor face 
бил one basis Of upstream distortion profiles alone. 
As discussed in Ref. 10 the upstream conditions are strongly 


dependent on the presence of the blade row itself. Carta and 
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Figure 19. Bent Vortex Model for Heat Addition. 
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Adamczyk present a graphic comparison of angle of incidence 
distribution as it differs immediately behind the distortion 
screen (upstream) from that at the blade row inlet. This 
angle of incidence distribution of Ref. 10 is compared in 
Fig. 20 with a semi-quantitative analysis of inlet flow 
ample based on the bent vortex sheet model for pressure dis- 
tortion. This analysis is discussed in Appendix B. 

There are no measured or predicted angle of incidence 
distributions available for comparison with the temperature 
distortion model. However, the encouraging agreement in the 
case of pressure distortion prompts a similar analysis of 
temperature distortion. Figure 19 shows the change in axial 
velocity at the points A, B, B', C, D and E, and Figure 21 
Shows the change in U, at these points. A semi-quantitative 
analysis similar to Appendix B gives the inlet flow angle 
(84) distribution shown in Figure 22. The area above the 
Hüminal clean inlet value of B1 Contribuie to stall. The 
large angles of 84 in Figurer? he ики илини 6” flow secs 
tor somewhat, but the majority of the high blade loading due 
te large values of B1 ika ии си 065 sector. Obsery— 
тмир ијә Illb), it can be seen that the compressor with 
dee corced flow stalls first, which gives qualitative verifi- 


Eatıon of the bent vortex model for temperature distortion. 
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VI. DISTORTION INDEX FORMULATION 


The results of Section V suggest an index to predict 
compressor instability or stall that combines total pressure 
wadal vorticity. Figure 21 indicates that the area 
experiencing the highest values of inlet flow angle are those 
Eu also wou NM experience posTt We TINOM vorticity. Since 
there are existing data and digital computer programs [Ref. 3] 


that compute radial vorticity as well as AP Chen t LS 


t? 
easy to search the compressor face for these values. There 
are several possibilities for processing these values, once 
they are calculated, into an index. One method would be to 
compare the sum of radial vorticity and AP, for a specified 
sector (percentage of compressor face area) with a critical 
value. Since there is no known critical value, another ap- 
proach that eventually might lead to a critical value 15 

recommended. In Figure 23 the total pressures across the 


face of the compressor are shown. P. is the average total 


pressure. If P. - P. 1s posıtıversriengehewerreet 1s 260 
Merca e he blade loading as shown in Figure 20. Positive 
radial vorticity also increases blade loading as shown in 
usure 203). and as in Refs. 1-3 and 6. The most signifi- 
cant effect on compressor instability or stall would be the 
Et value of AP, andain ad la vorticity. The pressure 
change can be normalized with Po, and an index can be formu- 
slated by summing, over a sector, the combined values of 2. 
and w 


po 
EN 








05 


2222 2222 222 
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Where the integer k represents the number of pressure rakes 
to be searched. By rotating or stepping the sector around 
the array of probes, it is possible to compute the maximum 
valued sector at any time. By plotting Ene maximum value 
of the index versus time, a graph similar to Figure 24 
could be analyzed for a possible critical value. A time 
history statistical comparison of deviation from the mean 
value of the index is another possible test to be applied 
and is done by Vann [Ref. 4] after adapting this proposed 
index to existing digital computer programs and comparing 
it with other known indices. 

It is implied that the index is a positive number at 
all times. There is nothing in the formulation that re- 
ти ити to be positive, and if the sector over 
which the index is computed is small enough then it may in 
fact be negative. Throughout this paper the term stall has 
been applied to large inlet flow angles that load the blade 
row beyond the stall limit of the individual blades. If 
the inlet flow angle is too small, a phenomenon known as 
negative stall occurs. Reference 8 discusses the phenomenon 
of nepative stall, but it is important to point out that a 
Irati endstortion index would not indicate a negative 


521011 the analysis of this index by Vann [Ref. 4], 
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Distortion Index 
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a sector of 60° was used, and the index values were always 
Positive. Reference 4 addresses the size of sector used to 


compute this index and other commonly used indices. 
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VII. CONCLUSIONS 


Interpretation of distorted inlet flow interaction with 
a compressor based on vorticity maps [Refs. 1-3] is satis- 
factory as long as the assumption of constant static pres- 
sure is valid. Experiment has shown that constant static 
Pr sure us not the case at the compressor face, and the 
ie model müst be expanded to account for this. The 
failure of the vorticity model to predict the conditions 
at the face of a compressor (and thus its loss of stall mar- 
gin) in the case of total temperature distortion is a result 
of not considering the effects of the compressor on the up- 
53051 107 113610. “he üpstream influence of the compressorT 
Causes an area change of the undistorted streamtube that re- 
sults in greater velocity changes than those induced by the 
ке исе The pressure ratio equation, Equation (37), 
and parallel compressor theory adequately predict trends in 
compressor stability. 

A model based on the bending of the vortex sheet gives 
results comparable to parallel compressor theory and in the 
case of total pressure distortion qualitatively duplicates 
Cartas results. 

MS cali incidence angle ratio suggests a distor- 
2511751 hit combines radial vorticity and total pressure 
gradients since both can be related to blade loadings that 


EX: ouch an index has been developed and has been 
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tested [Ref. 4] against other distortion indices using common 
data for comparison. The new distortion index has proven 
equal to or better than the distortion indices with which it 


was compared for accuracy in predicting compressor stall. 
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APPENDIX A 


Ingestion of hot gas from forward firing armament sys- 
tens has more than temperature effects on the performance 
of an axial-flow compressor. The characteristic line of 
corrected engine rotational speed (N//0 - constant) appears 
to be a function of “temperature alone since 0 is defined 
as the ratio of actual temperature to some reference tem- 
perature, usually a sea level standard day temperature. 
However, the reference to temperature alone is a convention 
that has evolved from the fact that compressors in gas 
turbine engines for aircraft applications use only air as a 
fluid medium. The corrected engine rotational speed (N//9), 
sometimes identified as referred speed, is a similarity 
parameter designed to ensure comparability of data taken 
under varying conditions of testing. In the case of a com- 
m0 blademrow it is desirable to preserve the blade 
angle of attack. Consider two "different" compressors as 
1071101010 €. ft the angle wer attack zs to remain the same 
for both compressors, then 8 must be equal to 8'. 8B is the 
inlet. flow angle but differs from the angle of attack only 


by a constant related to blade geometry as discussed in Ap- 





poudm sp E For similarity 
^ ! ! 
tan B = = -.. 
Z 7 
animas equal M', Therefore 
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17. 
aM a'M' 


where a is the speed of sound. Cancelling the Mach numbers 
and substituting N (usually RPM) for the angular speed Q, 


one obtains the relationship 


Ne Se 


a a 





meer Thezceompressor performance maps usually relate the 
ensrzcte ri Lie, 2ot a specıd Tevceoınpressar under differme 
conditions, the radius r could be eliminated and the similar 


rotational speeds expressed as 


1 
er U , or N=N' — 
a a a 
Since a^ = yRT, the ratio of sonic speeds reduces to the 


tatio of the square root of the temperature -if y and R are 
the same for both compressors. This is- usually the case 
Ance OTT for air within the temperature range 
encountered and y = y' is fundamental to similarity. The 
gas constant R varies with molecular weight, however; and 
if the fluid medium is air plus foreign gas (steam or pro- 
dicke 007 üs £ On for example), then the similarity rela- 
tionship must retain sonic velocity as a variable. In fact 
onc lisenmaitune. (in particular older literature) uses the 


emee ed FOtdtional speed notation of 


N 2 175775 
7: all 


65 





Ls awsHaracteristic length, usually the tip radius 
of the first stage if comparing different sized compressors. 
All corrected rotational speeds in this paper are the 
more usual N//Y@ since foreign gas distortion is not con- 
sidered. References 2 and 3 analyze the effects of foreign 
gas distortion or vorticity and how it affects compressor 


performance. 
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APPENDIT B 


Carta and Adamczyk [Ref. 10] define a. stall angle 
parameter c as the ratio of instantaneous incidence angle 
to the stall incidence angle of the compressor blade. They 
further define the incidence angle a to be the difference 
between the stagger angle A and the inlet flow angle By. 


From Figure 14 


a = By 27770t B1 - Qt A 


and B1 differs from a only by a constant ” 1s a func- 
tion of blade geometry. It will not be inconsistent then 
to plot a distribution of inlet flow angle 81 vs 8° instead 
ol o 520 as in Ret. 10. 

In this semi-quantitative analysis a nominal clean in- 
let (undistorted) value of 45° will be assumed for 51. The 
changes in inlet flow angle at the points A, B, B', C, D 


and E of Figure 17, using the change in U. and U, of Figures 


0 
172077 7 “rc predicted graphically ın Figure B-1. The 
Bigerineeecsinlet flow angles are plotted in Figure 21 with 
the stall parameter co from Ref. 10. 

Although the magnitude of many of the variables in this 
semi-quantitative analysis are somewhat arbitrary, the di- 
rections in which they were applied have been developed 


ки ene text of this paper. Deviations from the 


magnitude used should give the same general shape of inlet 
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Figure B-1. Compressor Effects on Inlet Flow Angle. 
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1707776 curve as that of Figure 21. The important point 

is that the general trend of inlet flow angle closely parallels 
that calculated by Carta and Adamczyk and reinforces the bent 
vortex model as having accounted for the upstream effect of 


the compressor on the flow field. 
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